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Abstract 

Using a collective coordinate numerical optimization procedure, we construct ground-state con- 
figurations of interacting particle systems in various space dimensions so that the scattering of 
radiation exactly matches a prescribed pattern for a set of wave vectors. We show that the con- 
structed ground states are, counterintuitively, disordered (i.e., possess no long-range order) in the 
infinite-volume limit. We focus on three classes of configurations with unique radiation scatter- 
ing characteristics: (i) "stealth" materials, which are transparent to incident radiation at certain 
wavelengths; (ii) "super-ideal" gases, which scatter radiation identically to that of an ensemble 
of ideal gas configurations for a selected set of wave vectors; and (iii) "equi-luminous" materials, 
which scatter radiation equally intensely for a selected set of wave vectors. We find that ground- 
state configurations have an increased tendency to contain clusters of particles as one increases the 
prescribed luminosity. Limitations and consequences of this procedure are detailed. 

PACS numbers: 
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I. INTRODUCTION 



A fundamental problem of statistical mechanics is the determination and understanding 
of classical ground states of many-particle systems - the zero-temperature particle arrange- 
ment that minimizes potential energy per particle. Although perfect crystalline (periodic) 
structures are often energetically favorable among all configurations, there is incomplete 
mathematical and intuitive understanding of the formation of order at low temperature,- 
introducing the possibility of disordered ground states. A disordered many-particle system 
is one that lacks long-range order. More precisely, disordered systems have a pair correlation 
function (72(1") (defined below) that decays to unity faster than |r|~'^~'^, for spatial dimen- 
sion d and some positive e, in the infinite-volume limit.- Recently, a collective coordinate 
approach has been used to identify certain possibly disordered ground states.-"^*^ Although 
the previous studies are suggestive of the relation between disordered ground states and col- 
lective coordinates for finite systems, a systematic investigation of disordered ground states, 
including whether they exist in the infinite-volume limit, is not yet available. 

In this paper, we use an "inverse" approach to construct classical disordered ground 
states with precisely tuned wave scattering characteristics via the aforementioned collective 
coordinate procedure. In a recent example of a "forward" problem, the scattering from glass 
ceramics with nanometer-sized crystals was likened to that of random sequential adsorption 
(RSA) of hard spheres,- a well-known, disordered many-body configuration.^ These ceramics 
are of interest in photonics applications because they are mechanically rigid and nearly 
suppress all scattering at long wavelengths.- 

By contrast, our method utilizes an inverse approach: we prescribe scattering character- 
istics {e.g., absolute transparency) and construct many-body configurations that give rise 
to these targeted characteristics. Potential applications include designing ground-state ma- 
terials as radiation filters or scatterers, and materials transparent to specific wavelengths of 
radiation, among others. We apply our methodology initially for structureless (i.e., point) 
particles. However, it could be generalized to structured particles, colloids, or as bodies using 
the appropriate structure factors for finite-sized particles as is done for random media.- 

Since previous studies utilized small periodic systems,-i^i^ we first establish that system- 
atically increasing the system size has no effect on the degree of disorder. Extrapolation 
from these results indicates that the constructed configurations remain disordered in the 
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infinite-volume limit; a seemingly counterintuitive proposition. We then construct disor- 
dered ground states with special scattering properties: "stealth materials," "super-ideal 
gases," and "equi-luminous materials." 

We use the term "stealth" materials to refer to many-particle configurations that com- 
pletely suppress scattering of incident radiation for a set of wave vectors, and thus, are 
transparent at these wavelengths.— Periodic (i.e., crystalline) configurations are, by defini- 
tion, "stealthy" since they suppress scattering for all wavelengths except those associated 
with Bragg scattering. However, we construct disordered stealth configurations that prevent 
scattering only at prescribed wavelengths with no restrictions on any other wavelengths. 

We define a "super-ideal" gas as a single many-particle configuration whose scattering 
exactly matches that of an ensemble of ideal gas configurations, or Poisson point distribu- 
tions, for a set of wave vectors. A super-ideal gas has a structure factor that is identically 
unity for specified wave vectors, and thus, this single configuration would be impossible to 
differentiate from an ensemble of ideal gas configurations for the specified wave vectors. 

We define an "equi-luminous" material to be a system whose scattering is constant for 
a set of wavelengths. Although stealth materials and super-ideal gases are subsets of equi- 
luminescent materials, we use this term to refer to materials that scatter radiation more 
intensely relative to an ideal gas. These materials that scatter radiation much more intensely 
than an ideal gas for a set of wave vectors have enhanced density fluctuations and show local 
clustering similar to polymers and aggregating colloids. Typically, scattering experiments 
on these systems are used to shed light on the characteristic length scales of a system.— >^>^ 
With our inverse procedure, we impose the degree of clustering by tuning the scattering 
characteristics for certain wavelengths. 

Upon generating ensembles of ground-state configurations for each class of materials 
described above, we characterize local order of each ensemble. We place emphasis on pair 
information in real space via the pair correlation function (72(1") and in reciprocal space 
through the structure factor S'(k) as these functions are experimentally accessible and used 
widely in many-body theories.-»^>^ 

The pair correlation function g2{r) is the normalized two-particle probability density 
function P2(r) and is proportional to the probability of observing a particle center at r 
relative to a particle at the origin.- For a statistically homogeneous and isotropic medium, 
the pair correlation function (72(1") depends only on the magnitude of r = |r|, and is commonly 
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referred to as the radial distribution function g2{'r), which henceforth is the designation used 
in this paper. 

The structure factor S'(k) is proportional to the intensity of scattering of incident radia- 
tion from a configuration of particles and is defined as 

S(k) = (1) 

where p(k) are the collective coordinates and k are the wave vectors associated with the sys- 
tem volume and boundary conditions. Collective coordinates p(k) are the Fourier coefficients 
in the expansion of the density field: 

N 

p(k) = ^exp(«k-r,), (2) 

j=i 

where denotes the location of particle j. When 5'(k) depends only on the the magnitude 
oi k = |k|, the structure factor S{k) is related to the Fourier transformation of g2{r) — 1, 
ignoring the forward scattering associated with k = 0, 

Sik) = '^ + pj exp(ik ■ r) [^2(r) - 1] dr, (3) 

where p is the number density. For highly ordered systems, both g2{r) and S{k) contain a 
series of 5-functions or peaks at large r and k, indicating strong correlations at the associated 
pair distance. In configurations without long-range order, both g2{r) and S{k) approach 
unity at large r and k. 

Several inverse methods have sought to construct systems using pair information 
in real space, particularly in addressing the question of pair correlation function 
"realizability,"-'^»^»^»^'^ which asks whether a given pair correlation function, at num- 
ber density p, can be realized by spatial arrangements of particles. Typically, in real space, 
a "target" radial distribution function is chosen and many-particle configurations are found 
that best match the target g2{f)- Stochastic optimization techniques have been a popular 
reconstruction method in finding spatial arrangements of particles that best approximate a 
target correlation function.— i^i2^i22i'2i 

In contrast to these real-space methods, we target pair information in reciprocal space 
to construct configurations whose structure factor exactly matches the candidate structure 
factor for a set of wavelengths. In addition, our procedure guarantees that the resulting 
configuration is a ground-state structure for a class of potential functions. 



The remainder of this paper is as follows. Section [TTl presents background on disordered 
ground states and motivates our choices for candidate structure factors, while Section [XTTl 
outlines the numerical procedure. Structure factors, radial distribution functions, and rep- 
resentative particle patterns for stealth materials, super-ideal gases, and equi-luminous ma- 
terials are found in Sees. HVl IVl and |Vll Lastly, general conclusions and discussion relevant 
to classical disordered ground states and this procedure are found in Sec. I VI II Appendix |X] 
compares minimization algorithms and analyzes the energy landscapes associated with our 
potential functions. 

II. STEALTH MATERIALS, SUPER-IDEAL GASES, AND EQUI-LUMINOUS 
MATERIALS 

Physical intuition and experimental facts suggest that in the zero-temperature limit, 
classical systems of interacting particles adopt a periodic structure to minimize potential 
energy. The "crystal problem" has attempted to determine the fundamental mechanism 
that forces particles into ordered states, but the existence of these mechanisms has yet 
to be fully understood. - The notion of disordered ground states is particularly mysterious 
because of the lack of symmetry, lack of long-range order, and degeneracy of ground-state 
configurations. 

The characterization of order in solid phases in the low-temperature limit has been well 
studied. In addressing the crystal problem, it has been suggested that nonanalyticity of 
thermodynamic functions may yield "turbulent," or nonperiodic, Gibbs states at positive 
temperature.-SS As a consequence, a turbulent crystal, characterized by fuzzy diffraction 
peaks, is possible as a nonperiodic solid phase, in addition to periodic and quasiperiodic 
structures. Turbulent crystals have been examined previously^^i^I and evidence has been 
presented that at low temperature, equilibrium states may contain disorder.— Theoretical 
work has created classical lattice models with short-range interactions whose ground states 
contain the property of disorder.— In addition, a simple gradient model was used to develop 
a disordered state, a labyrinth, within a pattern forming system. Although this state was 
ultimately excluded from being a ground state, the authors were unable to exclude other 
models as potentially yielding disordered ground states.^ Despite significant research at- 
tention, understanding of the attainability of disordered configurations as classical grounds 
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states is incomplete. 

We choose to limit this investigation to the study of stealth materials, super-ideal gases, 
and equi- luminous materials based on the potential applications and fundamental interest. 
Constructing systems with transparency at specified wavelengths, including wavelengths 
outside of the low-A; region, is the primary motivation. A stealth material has a structure 
factor that is exactly zero for some set of wavelengths. In the \ow-k limit, several disordered 
systems nearly suppress all scattering, though most are not ground states. For example, 
when crystallizing hard colloidal spheres, stacking faults of close-packed layers create de- 
viation from perfect crystallinity. The structure factor is nearly zero in the low-A; region 
followed by a strong Bragg peak and a diffuse peak before decaying to unity.— 

The examples of stealth configurations can also be described as "hyper uniform." Hyper- 
uniform systems have the property that 

lim S{k) = 0, (4) 

i.e., infinite- wavelength density fluctuations vanish.— Hyperuniform point patterns arise in 
the structure of the early Universe,— maximally random jammed packings,— certain 
tilings of space,— and the ground-state configurations of certain repulsively interacting par- 
ticle systems.— 

The structure factor for a super-ideal gas is exactly unity for a set of wave vectors and 
unconstrained for the remaining wave vectors. We narrow our study to constraining wave 
vectors in the small k region. We choose super-ideal gases based on interest in Poisson 
point distributions. The Poisson point process has the simplest candidate S{k) and g2{f), 
both being exactly unity for all k and r, respectively. For a single finite configuration, 
the structure factor exhibits random fluctuations about unity. But in a super-ideal gas, 
we constrain a set of wave vectors so that S{k) is exactly unity and, for all wave vectors 
outside of this unconstrained set, the intuitive expectation is that the structure factor would 
ensemble- average to unity. In fact, this does not necessarily happen and an interesting 
alternative behavior arises, discussed in detail below. 

"Equi-luminous" describes materials that scatter light equally intensely for a set of wave 
vectors. The structure factor for this class of materials is simply a constant for a set of wave 
vectors. Subsets of equi-luminous materials include super-ideal gases {S{k) = 1) and stealth 
materials {S{k) = 0). Here, we focus on materials whose structure factor is a constant 
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greater than unity, as these systems show strong local clustering and intense scattering 
relative to an ideal gas. 



III. NUMERICAL PROCEDURE 

The numerical optimization procedure follows that of Uche, Stillinger and Torquato^ 
used to tailor the small k behavior of the structure factor. The structure factor S{k) and 
collective coordinates p(k), defined in Eqs. ([T]) and ([2]), are related to the quantity C(k), 

5(k) = l + |c(k), (5) 

where 

N-l N 

^(k) = 5^ 5^ COS [k ■ (r, - r,)] . (6) 

j=i i=j+i 

For a system interacting via a pair potential f (rj — rj), the total potential energy can be 
written in terms of C(k), 

= n-'Y.^ik)cik), (8) 

k 

where Q is the system volume and V^(k) is the Fourier transform of the pair potential function 

V{k)= / rfrf(r)exp(ik-r). (9) 
Jn 

For a region of space with dimensions L^, x Ly, or Lr,. x Ly x in one, two, or three 
dimensions, subject to periodic boundary conditions, the infinite set of corresponding wave 
vectors has components 

= ^. (10) 

where are positive or negative integers, or zero and 'j=x, y, z as needed. For example, in 
three dimensions, the set of wave vectors are 



27™^; 27™,„ 271712 



Lx Ly Lz 



(11) 



We introduce a square mound V^(k) that is a positive constant Vq for all k G Q, where 
Q is the set of wave vectors such that < |k| < K, and zero for all other k. In the 
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infinite- volume limit, this corresponds to a system of particles interacting via a real- space 
pair potential function that is bounded, damped, and oscillating about zero at large r.->^ For 
a cutoff radius K, there are 2M{K) wave vectors in the set Q, where M{K) is the number 
of independently constrained collective coordinates. That is, constraining C(k) implicitly 
constrains C(— k) due to the relation 

Cik) = C{-k). (12) 

For a system of particles in d dimensions, there are dN total degrees of freedom. We 
introduce the dimensionless parameter x to conveniently represent the ratio of the number 
of constrained degrees of freedom relative to the total number of degrees of freedom 

^ = -ir- 

The global minimum of the potential energy defined in Eq. ([8]) has the value of 



^ ' kGQ 

if and only if there exist particle configurations that satisfy all of the imposed constraints, 
which necessarily occurs for x < 1. Minimizing Eq. ([8]) to its global minimum, for % < 1, 
yields ground-state configurations that are stealthy for all kGQ. 

To target a specific form of the structure factor to certain nonzero values, such as S{k) 
= 1, we introduce a second nonnegative objective function, 

$ = 5^V(k)[C(k)-Co(k)]% (15) 

where Co(k) is associated with the target structure factor by Eq. ([5]). If Eq. (1151) is taken to 
be the potential energy of an A^-body system, then two-, three-, and four- body interactions 
are present.- Equation ( fTSl) has a global minimum of zero, for x < 1, if and only if there 
exist configurations that satisfy all of the imposed constraints. Minimizing Eq. ( |T5i) is used 
to construct super-ideal gases and equi-luminous materials as ground-state configurations. 

Three algorithms have been employed previously for minimizing Eqs. ([8]) and (fT5l) : steep- 
est descent,- conjugate gradient,^ and MINOP.^ Steepest descent and conjugate gradient 
methods are line search methods that differ only in their choice of search directions.— The 
MINOP algorithm is a trust-region method. When far from the solution, the program 
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chooses a gradient direction, but when close to the solution, it chooses a quasi-Newton 
direction.— Upon each iteration, the program makes an appropriate update to approxi- 
mate the Hessian.— 

We find that that neither the conjugate gradient method nor MINOP algorithm signif- 
icantly biases any subset of ground-state configurations. The resulting configurations are 
visually similar, and the ensemble-averaged radial distribution function and structure factor 
produced by both methods have similar features. We chose the MINOP algorithm because 
it has been demonstrated to be better suited to the collective coordinate procedure than the 
conjugate gradient method.- We refer the reader to Appendix lAl for characterization of the 
energy landscape and comparison between line search methods and MINOP. 

Three sets of initial conditions were considered: random placement of particles (Poisson 
distributions), random sequential addition (RSA), and perturbed lattices (integer, triangu- 
lar, and face centered cubic in one, two, and three dimensions respectively). For an RSA 
process, particles are assigned a diameter and randomly and irreversibly placed in space such 
that particles are not overlapping.^ At sufficiently high x, usually x ^ 0-6, the constructed 
ground-state systems apparently lose all memory of their initial configurations. The analyses 
presented in the following sections will be those of random initial conditions. In some cases 
at large x, a global minimum is not found. For the results discussed here, Eqs. ([8]) and (fT5|) 
were minimized to within 10~^^ of their respective minimum value. All other trials were 
excluded from the analysis. 

The region of space occupied by the N particles was limited to a line in one dimension, a 
square in two dimensions, and a cube in three dimensions, with periodic boundary conditions. 
For stealth materials, particular attention was paid to the choice of N for two and three 
dimensions. Minimizing Eq. ([8]) for large x is known to yield crystalline ground states.-i^ 
We choose to be consistent with previous studies. In two dimensions, was chosen as a 
product of the integers 2pg, and p/g is a rational approximation to 3^/^ so that all particles 
could be placed in a triangular lattice configuration without substantial deformation. In 
three dimensions, was usually chosen so that A^ = 4s'^, where s is an integer, so that the 
particles could be placed in a face centered cubic lattice without deformation. In minimizing 
Eq. (fT5|) . A^ occasionally was assigned other values. 
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IV. RESULTS FOR STEALTH MATERIALS 



A. Infinite- Volume Limit 

Previous work suggesting the existence of disordered ground states utilized small simu- 
lation boxes containing up to several hundred particles.-'^'^ Our goal here is to show that 
constructed systems continue to show no long-range order in the infinite- volume limit. For 
d=2 and 3, systems containing up to several thousand particles were constructed by mini- 
mizing Eq. ([8]) for small x values. We find that ground-state configurations differing only 
by A^, with p and x fixed, are disordered and exhibit the same local structure. 

Figures [T] and [2] demonstrate the behavior of S{k) and g2{T) for stealth materials con- 
strained at x=0.05. The structure factors for systems only differing in N have identical 
characteristics. The structure factor is exactly zero for constrained wave vectors and subse- 
quently peaks above and fluctuates about unity. The averaged S{k) initially peaks to a value 
of 1.10 and decays rapidly to unity, a feature that is more apparent for a large ensemble of 
conflgurations. 

The radial distribution function remains essentially invariant as the number of particles in 
the simulation box increases from 108 to 6912, with x and p flxed, as shown in Figure [2l For 
a single realization, a system containing 108 particles has the same shape of g2{'r') as that of 
the larger system but shows signiflcant statistical noise. In the flgure, we ensemble-average 
the results for the smaller system to make clear the structural similarities. In the smaller 
system, the large-r behavior is unavailable due to the minimum image convention of the 
periodic box. Thus, the flgure only displays local structure, which is clearly disordered. In 
both cases, g2{r) dips slightly below unity for small r and quickly approaches and oscillates 
about unity with a diminishing amplitude. 

B. Effect of Increasing Constraints 

It was previously reported that minimizing collective density variables for sufficiently 
high X induces crystallization,-'^'^ therefore, for disordered, stealth ground states, we have 
minimized Eq. ([8]), focusing on the low x regime. 

For (i = 1, crystallization occurs for x > 0.5.^ For most values of x below the crystalliza- 
tion threshold, the structure factor is zero for all constrained wave vectors, it then peaks 
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above unity immediately outside K and decays toward one. As x is increased, the height of 
the peak decreases and at x higher than 0.30, a second peak forms. Below the crystallization 
value, S{k) dampens to unity for k larger than K. 

For d = 2, three x regimes have been reported: disordered for x <0.57, "wavy crystalline" 
for 0.57 < X < 0.77, and crystalline for x > 0.77 A We choose to investigate well below these 
ordered regions. In the ensemble-averaged structure factor, a peak forms in S{k) immediately 
beyond K and decays toward unity. As x increases, the magnitude of the peak increases, in 
contrast to d = 1, and at the wavy crystalline threshold, several peaks begin to form beyond 
K. The height of the peak is density dependent, however, S{k) generally has a maximum 
between 1.1 — 2.0 in the disordered region. 

For d = 3, the transition from disordered to crystalline regimes was identified previously 
to occur at x near 0.5.^ Constraining x below 0.45, S{k) peaks immediately beyond K 
and smoothly decays to unity, but at x = 0.45467, S{k) smoothly oscillates about unity. 
The magnitude of the peak is generally smaller than for systems similarly constrained in 
lower dimensions. In Figure [3l we compare the structure factor for several x values for 500 
particles in a unit cube. We include a nearly crystalline system of 500 particles constrained 
at x=0. 54867. The order of the system is apparent by the series of sharp peaks in S{k) that 
persist at large k. 

We find that particles have a repellent core that increases in strength with increasing x- 
Figure m demonstrates the repellent core effect via the radial distribution function associated 
with the 500 particle system described above. For this particular system, at x = 0.45467, 
an exclusion region develops where g2{r) is exactly zero for a region near the origin. At x 
= 0.54867, the peaks demonstrate crystallinity. 

Increasing x tends to increase the the net repulsion of the potential, which is clearly 
observed in particle patterns. Since differentiating between disordered stealth systems is 
most instructive in two dimensions, we present particle patterns in this dimension only. 
Figure [5] compares particle patterns of 168 particles with Eq. ([8]) constrained for small x- The 
circular window in the figure represents the length scale of the wavelength associated with 
K. At the lowest x considered, the particles do not appear to have any spatial correlation. 
At higher x values, particles develop an exclusion shell about their center but do not have 
any long-range order. 
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C. Stealth Materials Spherically Constrained in k Space 

The entirety of disordered stealth materials studied have involved constraining wave vec- 
tors in a spherical shell near the origin. Since this procedure is capable of constraining wave 
vectors of choice, we have constructed configurations in which two disconnected, concentric 
regions of wave vectors are constrained, one near the origin and a second further from the 
origin. 

We define parameters Kq < Ki < K2 < K3 as magnitudes of limits for constrained 
wave vectors. In this class of stealth configurations, we constrain collective coordinates so 
that S{k) is zero for all wave vectors in two spherical shells about the origin in reciprocal 
space. Specifically, S{k) is zero for all Kq < |k| < Ki and all K2 < |k| < K^,. The region 
-^1 < |k| ^ -^2 is defined as the intermediate region, where the structure factor can be free 
to fluctuate or to be controlled. Figure [6] shows the location of parameters. These systems 
are constructed by introducing a l^(k) that contains two square mounds. 

Spherically constrained stealth configurations were constructed in two and three dimen- 
sions. We present ensemble-averaged radial distribution functions and structure factors for 
d = 3 and particle patterns ior d = 2 as these provide the clearest representation of the 
general trends. 

The peaking phenomenon in S{k), as originally observed in simple stealth materials, is 
evident both immediately beyond Ki and beyond K^. The structure factor increases above 
unity slightly in the intermediate region but peaks and decays rapidly to unity beyond K3. 
This is seen for all test cases and it is not immediately clear if this phenomenon persists 
for large separations of constrained regions (i.e., K2 — Ki). Figure [H] shows the ensemble- 
averaged radial distribution function and structure factor for configurations of 500 particles 
constrained so that S{k) = for < |k| < S.Stt and IStt < |k| < 14.87r. 

An important feature of this procedure is the ability to suppress scattering for wave 
vectors that are normally Bragg peaks in a crystalline material. With a 500 particle, three- 
dimensional system, minimizing Eq. (JHl) for a single square mound V(k) for x = 0.54867 
creates a crystalline ground state where the first Bragg peak occurs just beyond k = lAn. 
The constructed system in Figure [6] suppresses scattering for a range of k surrounding 147r. 
Because the intermediate set is free to fluctuate, the total number of constrained wave vectors 
in this stealth configuration is less than that of a crystalline configuration. The stealth region 
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away from the origin can be shifted to larger k and can be tailored in magnitude. 

The scattering of the intermediate set of wave vectors can also be controlled using Eq. (|T5l) . 
One such configuration that can be developed is a stealth/ideal gas hybrid. In these many- 
particle systems, there is behavior typical of crystalline solids at selected regions of reciprocal 
space and ideal gas at others. In Figure [3, a 500 particle system in three dimensions is 
constrained so that S{k) = for < |k| < S.Svr and ISvr < |k| < 14.87r. The intermediate 
set is constrained to S{k) = 1. 

The characteristics of the radial distribution function vary depending on the constraints 
of the intermediate set. In all cases, the radial distribution function shows weak neighbor 
peaks before approaching unity. The second neighbor peak, though, is stronger than the 
first neighbor peak, a trait uncommon to most conventional many-body systems. The extent 
of the repelling-core region varies depending on the chosen value of Ki and the extent to 
which the intermediate set is controlled, evidenced by Figures [6] and [71 Further reducing the 
the value to which the intermediate set is controlled is likely to increase the repelling-core 
region given the relation between S{k) and g2{r) in Eq. ([3]). 

Spherically constrained stealth particle patterns consisting of 168 particles with the in- 
termediate set respectively unconstrained and constrained are shown in Fig. [HI To serve as a 
basis of comparison, a realization of a wavy crystalline configuration generated by suppress- 
ing scattering for all wave vectors up to 227r is also shown. The wavy crystalline material, 
Fig. [HH, stands in sharp contrast to the stealth materials since particles tend to align in 
well-patterned strings. The spherically constrained stealth materials lack any order. With 
the intermediate region uncontrolled. Fig. [8b, particles tend to align into weak "strings." 
Creating a stealth/ideal gas hybrid. Fig. [Ht, decreases the tendency to align in weak strings. 
The diameter of the circular window is equivalent to the length scale of K3, indicating the 
ability to impose system features with a specified length scale. 

V. RESULTS FOR SUPER-IDEAL GASES 

Super-ideal gases were constructed at various x values for (i = 1, 2, and 3. The maximum 
attainable x value varied depending on spatial dimension and system size but was generally 
near x = 0.95 for most initial configurations. At x near unity, the minimization routine 
sometimes failed to find a global minimum of $. However, in ci = 3 and x near unity, the 
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success rate for finding global minimum was much improved over that of lower dimensions. 

The results from all dimensions studied have similar characteristics. For a single realiza- 
tion, the structure factor is exactly unity for k < K. Outside of the constrained region, the 
structure factor seemingly fluctuates about unity. However, for an ensemble of super-ideal 
gases, the structure factor had a small peak immediately beyond K that slowly decayed to 
unity. This small peak is unexpected since we impose no constraints on S{k) for wave vec- 
tors outside K and would expect that correlations not exist at these k. The construction of 
super-ideal gases reveals a subtle coupling between S{k) within the constrained region and 
S{k) outside the constrained region that manifests upon ensemble averaging. For d = 1, 
2, and 3, S{k) never exceeds a value of 1.25, 1.18, and 1.10 respectively when ensemble 
averaged. For x < 0-4 and x > 0.96, the peak generally decays to unity rapidly. For 
0.4 < X ^ 0.96, the peak is rather long ranged, decaying much more slowly. Figure [9] shows 
the ensemble- averaged g2{'r') and S{k) for two systems containing 500 particles in = 3. 

The radial distribution function has characteristics common across dimensions studied. 
For r > O.lLj:, (?2(r) shows very small fluctuations about unity. The local behavior of g2{r), 
r < O.lLx, varies depending on x and is most sensitive to x in = 3. Figure [9] shows that 
super-ideal gases at x = 0.90667 exhibit severe local clustering as the radial distribution 
function has a contact value (72(0) near 7. However, at x = 0.98967, the structure more 
closely resembles an ideal gas. Figure [TOl tracks the contact value of g2{r) for various x- At 
significantly large x, local clustering is suppressed and the super-ideal gas structure closely 
resembles that of an ensemble of ideal gas configurations. An interesting consequence is that 
the local structure of a super-ideal gas at very small x resembles that of a super-ideal gas 
at very large x- 

Two-dimensional particle patterns reveal subtle differences between super-ideal gas con- 
figurations and Poisson distributions. Figure [11] compares a realization of Poisson distri- 
bution of 418 particles and a super-ideal gas at x=0.90. At certain x, super-ideal gases 
exhibit local order and a tendency to align into weak "strings" that is best revealed only 
upon ensemble averaging. At smaller x, super-ideal gases are particularly difficult to discern 
from a Poisson distribution. 
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VI. RESULTS FOR EQUI-LUMINOUS MATERIALS 



Equi-luminous materials represent a broad class of materials that scatter light equally 
intensely for a set of wave vectors, where a super-ideal gas is a special case. We choose to 
focus here on equi-luminous materials whose scattering in the small- /c region is much more 
intense relative to that of an ideal gas. For all cases considered, we observe qualitative 
similarities in the ensemble-averaged radial distribution function and structure factor. The 
structure factor is exactly equal to the chosen constant for all k < K. Beyond K, the 
ensemble- averaged structure factor decays to unity very slowly at a rate dependent on the 
specified constant in constrained region. Figure [T^ compares S{k) for several equi-luminous 
configurations containing 168 particles for d = 2. These configurations have constrained 
S{k) values that are exactly 1,2,3, and 4, respectively, for x = 0.34523. Achieving x above 
0.37 for S{k) = 4.0 for = 2 proved challenging as the minimization procedure often failed 
to find global minima. 

In real space, we observe that constraining S{k) to be increasingly large has no affect on 
the long-range behavior of (72 (^)- Generally, for r > O.IL^, g2{r) averages near unity. How- 
ever, for r < O.IL^, strong local correlations rapidly vanish for increasing r. By increasing 
the constrained value of S{k), the contact value of g2{r) increases significantly but does not 
change the large r behavior, which remains near unity. Figure [T3] demonstrates the behavior 
of g2{r) corresponding to the systems described above. 

Realizations of equi-luminous materials demonstration the aggregation of particles com- 
mon to this class of equi-luminous materials. Figure shows configurations for which S{k) 
= 2 and S{k) = 4 for x = 0.37. Particles tend to cluster and align into well-formed strings, 
or filamentary structures, as opposed to clustering radially. Filamentary structures arise 
in astrophysical systems, particularly for the distributions of galaxies. For a larger target 
S{k), the aggregation is increasingly severe and particles nearly stack on top of each other. 

In the extreme case of targeting S{k) toward its maximum value of for very small 
X, particles tend to collapse upon each other, yielding an overall rescaling of the system 
length scale. It should be noted that in these extreme cases, global minima were rarely 
found and required successive iterations to achieve optimality. For example, in seeking a 
ground-state configuration with S{k) = 10 for x = 0.01785, the ground-state configuration 
in which S{k) = 8 for x = 0.01785 was used as the initial condition. 
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VII. CONCLUSION AND DISCUSSIONS 



Previous studies concerning constraints on the collective coordinates of particle 
configurationa^i^'^ were suggestive that potentials defined by Eqs. ([H]) and ffT^ yield classical 
disordered ground-state configurations. These studies restricted the system size to small 
periodic boxes consisting of at most 12 particles in d = 1/ 418 particles in d = 2,^ and 500 
particles in ci = 3.- In this investigation, we find that increasing the system size while fixing 
p and X does not affect the structure of resulting configurations, suggesting that there are 
no long-range correlations in the infinite- volume limit. Constructing ground state materi- 
als via MINOP for systems significantly larger than that presented in Section IIVI becomes 
computationally challenging and currently is the limitation on system size. Regardless, our 
numerical evidence for disordered ground states in the infinite-volume limit corroborates 
recent analytical work that suggests the existence of energetically degenerate and aperiodic 
ground states in the infinite- volume limit via a similar potential.— 

Three novel classes of ground-state materials with potential radiation scattering applica- 
tions have been introduced: stealth materials, super-ideal gases, and equi-luminous materi- 
als. Each provides an unique opportunity to impose an underlying structure for a ground- 
state configuration with a known potential. With stealth materials, we have the precise 
control to suppress scattering at specified wavelengths. For a single square-mound ^(k), 
increasing the number of constraints on wave vectors near the origin drives systems toward 
crystallinity.^''^'^ However, by introducing V(k) with two square mounds and choosing to 
suppress scattering in two disconnected regions of reciprocal space, we disrupt the drive 
toward crystallinity. Particle patterns then have strong local correlations with an imposed 
length scale, but lack long-range order. It is interesting to note that all stealth ground-state 
configurations generated in this investigation are also hyperuniform point patterns; i.e., con- 
figurations that suppress density fiuctuations in the infinite- wavelength limit. In choosing 
which regions of reciprocal space to be stealthy, we exhibit the unique ability to control and 
suppress density fiuctuations at specified wavelengths. 

Equi-luminous materials investigated in Section |Vl] demonstrate ground-state structures 
that scatter light for a set of wavelengths more strongly than an ideal gas. In choosing 
the the number of constrained wave vectors and the value of S{k), one can design ground- 
state configurations that have increasingly clustered behavior for increasing x and S{k). By 
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adjusting we can impose correlations over a certain length scale. 

The differences between super-ideal gases and ideal gases are not intuitively obvious. 
The resulting local structure deviates significantly from that of an ideal gas as there is some 
degree of local clustering, with the exception of x near and x near 1. The propensity 
for super-ideal gases to exhibit local clustering suggests that, among the many energetically 
degenerate global minima in the landscape, configurations with some degree of clustering 
dominate relative to those that do not. We also find that this local order is particularly 
sensitive to the number of constraints imposed suggesting that the energy landscape changes 
significantly for small changes in x- Despite the local clustering, the resulting configurations 
exhibit no long-range order. 

One potential application of this procedure is in addressing the realizability question. 
Despite receiving significant attention, only necessary conditions have been placed on the 
pair correlation function (72 (r) and its corresponding structure factor S{k) for the realizability 
of a point pro cess General sufficient conditions have yet to be developed. This 

method can potentially address realizability from reciprocal space, which stands in contrast 
to real-space numerical reconstruction techniques.— i^ii^SiiSSii^ A limitation in its use for the 
question of realizability is the peaking phenomenon in S{k). In all cases studied, S{k) 
has a peak immediately beyond K, suggesting collective coordinates are not necessarily 
independent of each other and that the nonlinearity of the coordinate transformation plays an 
important role in constraining collective coordinates. Additionally, our results show that the 
magnitude of the peak is reduced for higher dimensions, suggesting that dimensionality plays 
a key role in the coupling among collective coordinates. The understanding of constraints on 
collective coordinates is incomplete. New questions arise: how are the wave vectors beyond 
K influenced by constraints below K and what role does dimensionality play in the peaking 
phenomenon? 

Another possible application of this procedure is to produce ground-state configurations 
that are ordered or quasiperiodic over a specified length scale for the design of photonic 
materials. Materials with a photonic band gap are of significant interest due to their tech- 
nological applications,—'^'^ and materials with desired band gaps have been designed.— 1^ 
From first-order perturbation theory, scattering of radiation is related to the photonic band 
gap of a material.-^ By targeting a structure factor that is maximal for certain wave vectors 
and zero for others, we may be able to construct ground-state structures with a specified 
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desired photonic band gap to a first-order approximation. 

In the present investigation, the potential ^(k) was chosen to be a square mound with 
compact support at K. However, the constructed ground-state configurations are equivalent 
to the ground-state configurations associated with a broad class of V^(k). If a ground-state 
configuration of density p is constructed by minimizing Eq. ^ or Eq. (IT51) to its global 
minimum value, Eq. (fT4l) or zero respectively, then the ground-state configuration is also 
a ground state at density p for any positive ^(k) that is bounded with the same compact 
support at K. 

Recently derived duality relations allow us to identify further degeneracies among ground 
states constructed via Eq. ([H]) and provide bounds on the minimum energy for various forms 
of f (r). These duality relations link the energy of configurations for a bounded and integrable 
real-space pair potential function to that of its Fourier transform and provide a fundamental 
connection between the ground states of short-range pair potentials to those with long-range 
pair potentials. More specifically, if a Bravais lattice structure is the ground state energy 
Umin for a pair potential function f (r), then Umin provides an upper bound on the ground 
state energy Umin of a system interacting via the Fourier-transform of the pair potential 
function V{\i) at the associated reciprocal density.— In using Eq. ([8]), we used a long- 
ranged, damped, and oscillatory t'(r), whose Fourier transform is a square mound l^(k). For 
X < 1 and at a given density p, this f (r) had a Bravais lattice as a ground-state structure, 
and for sufficiently large Xi ^ Bravais lattice becomes its unique ground-state structure at 
a certain density. By these duality relations, we know that the minimum energy of this 
Bravais lattice system, Eq. dHl), is the upper bound on the ground-state energy of systems 
at density p* = p^^{2n)~'^ interacting via ^(k).— 

Targeting specific forms of the structure factor with Eq. f|T5|) requires up to four-body 
interactions in real space to achieve these as ground states. However, it would be par- 
ticularly useful to construct disordered ground states using short-range, two-body poten- 
tials. Developing effective short-range, pair interactions, potentially via the Ornstein-Zernike 
formalism,— to achieve desired scattering properties remains a potential future direction. 
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APPENDIX A: ENERGY LANDSCAPE ANALYSIS 

For large x? the procedure sometimes terminates at a potential higher than the absolute 
minimum, failing to meet our criterion for global minimum. Generally this has been at- 
tributed to local minima in the energy landscape associated with Eqs. ([8]) and (ITSl) and the 
performance of the minimization algorithm.- Since little is known about these landscapes, 
we shed some light unto them and justify the discrepancy between the performance of the 
various algorithms. 

1. Energy Landscape 

The landscape of $ can be visualized for a system of three particles on a unit line. 
Imposing the constraint so that 5'(k)|k=-|-27r = j^D, the energy, $, becomes 

3 

$ = {cos[27r(xi - X2)] + cos[27r(xi - X3)] + cos[27r(x2 - X3)] + - - D}^. (Al) 

D must be within the range of realizability for S{k) (i.e., < D < The energy 

landscape possesses translational freedom, so we can fix xi at the origin, xi=0. Plotting Eq. 
flAll) versus particle coordinates for the case of D=0 and D=l provides a simple picture of 
the energy landscape. Figures [15] and [16] show the energy landscape and contour plot for 
D=0 and D=l respectively. 

For a stealth configuration, D = 0, Figure [T51 there is exactly one solution, corresponding 
to the crystalline arrangement of particles, that is a global minimum of $. The global 
maximum is a stacking of all particles onto the origin. For D > 0, Figure [16], there is a 
family of degenerate ground-state solutions located in a ring around configurational points 
corresponding to the periodic arrangement of particles. For this system, a single ground 
state exists only for D = 01 D = For all other Z), there exist degenerate ground 
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states. Taylor expansion about the minima for when D=0 indicates the landscape is quartic 
in nature. When D ^ 0, local maxima are located at configurational points associated 
with the periodic arrangement of particles, as evident in Figure Uni As D is increased, the 
local maxima are found increasingly higher in the landscape. Above D=1.5, the periodic 
arrangements of particle become global maxima and the origin becomes a local maximum. 

Visualization confirms that the energy landscape is nontrivial, particularly for D ^ 0. 
Local minima were not found despite our experience with the minimization procedure. Local 
minima can be found by constraining the next wave vector. For example, constructing a 
system for D = 1 at x=0.667 for a three-particle system introduces local minima. If 7^ 
and D 7^ $ cannot be minimized to zero at x=0.667. The crystalline arrangement 
achieves $=0 for D = and the stacking of particles achieves $=0 ioi D = 

2. Stationary Point Characterization and Capture Basins 

To enumerate and characterize the stationary points (local/global minimum or maximum, 
or saddle point) of the energy landscape $, Eq. ([8]), we introduce a nonnegative function, 
/(x), where x refers to particle coordinates: 

/(x) = V$ ■ V$. (A2) 

The global minima of /(x) have a value of zero and correspond to stationary points of the 
energy landscape $. Capture basins are the regions of phase space that, upon minimizing 
/(x), yield each stationary point. For three particles on a line, the stationary points are 
countable from the visualization of the landscape, and this provides a basis for comparison. 
We use steepest descent since it will generate nonoverlapping capture basins and we compare 
against MINOP. If at a certain point in the energy landscape MINOP's trust region overlaps 
with a different capture basin, it may favor a search direction that enters the overlapping 
capture basin and lead to a different stationary point. 

Table [T] summarizes the results of the capture basin analysis for a three particle one- 
dimensional system constrained at k=l and k=2 from several hundred trials. Random 
positions are assigned and /(x) is minimized, yielding a stationary point. Many times /(x) 
was not minimized to a global minimum indicating that a stationary point of $ was not 
found. In the table, the number of capture basins refers to the countable number of basins 
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in the landscape with one basin associated with each stationary point, without accounting 
for particle permutations. The fractions in the chart represent the fraction of the total 
number of stationary points found corresponding to each stationary point, or alternatively, 
the fraction of phase space that corresponds to a specific capture basin for each stationary 
point. 

The capture basins for a two-dimensional system of 12 particles were found for all avail- 
able values of Xi for D = 0. Because of the increased speed of MINOP, we only employ 
this algorithm to minimize /(x). Upon generating a stationary point, the potential $ was 
calculated and the eigenvalues of its associated Hessian matrix were found. Global minima 
of $ have nonnegative eigenvalues and $ = 0. Saddle points have at least one negative 
eigenvalue and $ > 0. Interestingly no local minima, local maxima, or global maxima were 
found. Table [TTl summarizes the capture basins and stationary points for 12 particles in 
two dimensions. Trials that did not yield stationary points are excluded from the reported 
fractions. 

No local maxima, global maxima, or local minima were found, however, they must exist 
in systems of sufficiently high dimension. The most unstable saddles generally have a high 
potential exceeding the global minimum by > 10^, while the most stable saddles have poten- 
tials $ that differed from a global minimum by 0(10^^). Our experience with finding local 
minima is consistent with these results. When MINOP terminates at $ above its global 
minimum value, which usually occurs when ^ 12, the potential of the final configuration 
is about 10~^ to 10~^ above the global minimum value. The strong correlation between the 
number of negative eigenvalues and potential is consistent with our experience with these 
algorithms. 

The number of saddle points grows rapidly with the total dimensionality of the system 
(number of particles times spatial dimension), and grows more rapidly than any other sta- 
tionary points. In the three-particle system at x=0.667, we observe that local maxima are 
located near saddle points in the energy landscape. It is likely that minimizing /(x) intro- 
duces a bias toward finding saddle points of $ rather than local maxima of Because local 
maxima are outnumbered by saddle points and the search procedure may favor directions 
toward saddle points of $, it likely that minimizing Eq. (]A2p would rarely generate local 
maxima for systems of high dimensionality. Global maxima are attained only by stacking 
particles in a single location. It is expected that this stationary point would be so greatly 
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outnumbered by all other stationary points that minimizing f{x) may never generate a 
global maximum. 

These simple studies allow us to better justify the inferences regarding the landscape 
associated with $. Overconstraining the three particle system for targeting a super-ideal 
gas introduced local minima. Increases to the system's dimensionality and to x increased 
the number of saddle points and is expected to increase the number of local minima. Com- 
parisons between MINOP and line search techniques indicate that the trust region often 
overlaps other capture regions which find global minima of $ more effectively. 
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TABLE I: Fraction of phase space corresponding to capture basins of landscape $ for three particles 
on a line as found by steepest descent (SD) and MINOP algorithms. 
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TABLE II: Fraction of phase space corresponding to capture basins of landscape $ for 12 particles 
in a unit square as found by the MINOP algorithm. 
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FIG. 1: Structure factor for stealth ground states for d=3, p=W8, and x=0-05. Increasing the 
system size from A^=108 to A^=6912 does not affect the scattering characteristics. The potential 
energy was minimized to within 10~^^ of its absolute minimum. 

FIG. 2: Radial distribution function for stealth ground states for d=3, p=108, and x=0.05. In- 
creasing the system size from Ar=108 to iV=6912 does not affect the resulting local structure. The 
potential energy was minimized to within 10~^^ of its absolute minimum. 

FIG. 3: Ensemble-averaged S{k) for stealth ground states consisting of 500 particles in a unit cube. 
At X > 0.45, S{k) begins to oscillate while damping to unity, (a) x = 0.11333, 250 reaHzations, 
(b) X = 0.25000, 50 realizations, (c) x = 0.45467, 50 realizations, (d) x = 0.54867, 4 realizations. 
The potential energy was minimized to within 10~^^ of its global minimum. 

FIG. 4: Ensemble-averaged g2{r) for stealth ground states of 500 particles in a unit cube. At x > 
0.45, g2{r) oscillates about unity with a shorter wavelength than observed for smaller x values, (a) 
X = 0.11333, 250 realizations, (b) x = 0.25000, 50 realizations, (c) x = 0.45467, 50 realizations, 
(d) X = 0.54867, 4 realizations. The potential energy was minimized to within 10^^^ of its global 
minimum. 

FIG. 5: Stealth particle patterns of 168 particles in two dimensions, (a) x = 0.04167, (b) x = 
0.20238. The bar below each graph and the circular window represent the characteristic length 
associated with K. Both systems are disordered but at higher x, particles tend to spread away 
from each other. The potential energy was minimized to within 10~^^ of its global minimum. 

FIG. 6: (a) Ensemble-averaged g2{r) and (b) S{k) for a stealth material of 500 particles in a unit 
cube. 25 reaHzations. S{k) = for < |k| < 8.87r and 137r < |k| < 14.87r and with the intermediate 
set unconstrained. The potential energy was minimized to within 10~^^ of its global minimum. 

FIG. 7: (a) Ensemble-averaged g2{r) and (b) S{k) for a stealth material of 500 particles in a unit 
cube. S{k) = for < |k| < S.Stt and 137r < |k| < 14.87r and with the intermediate set constrained 
to S'(fc) = 1. 10 realizations. The potential energy was minimized to within 10""^^ of its global 
minimum. 
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FIG. 8: (a) Wavy crystalline configuration generated by constraining S{k) = for all k < 22tt. (b) 
Stealth material generated by constraining S{k) = for < |k| < lOvr and 207r < |k| < 267r and 
the intermediate set is unconstrained, (c) Stealth material generated by constraining S{k) = for 
< |k| < IOtt and IStt < |k| < 227r and with the intermediate set constrained to S{k) = 1. The 
line beneath the figure and the circle in the figure approximately represent the length scale of K3. 
The potential energy was minimized to within 10~^^ of its global minimum. 

FIG. 9: Ensemble-averaged g2{r) and S{k) for a super-ideal gases in three dimensions. 500 particles 
in unit cube. % = 0.90667, 30 realizations, and % = 0.98933, 12 realizations. The dashed line shows 
location of unity for each structure factor. The potential energy was minimized to within 10~^^ of 
its global minimum. 

FIG. 10: The value of radial distribution function in the first bin for 500 particles in a unit cube. 
The contact value 52(0) increases initially for small x- However, for large X) local clustering is 
suppressed. Each data point represents at least 30 realizations, except for x = 0.98967 which 
represents 12 realizations. 

FIG. 11: Particle configurations of 418 particles, (a) Poisson point process, (b) super-ideal gas, 
S{k) = 1, X = 0.90. Ensembles of super-ideal gases reveal the presence of local clustering. 

FIG. 12: Ensemble-averaged S{k) for equi-luminous ground states consisting of 168 particles for 
d = 2. X = 0.34523, 50 realizations. The potential energy was minimized to within 10~^^ of its 
global minimum. 

FIG. 13: Ensemble-averaged g2{r) for super ideal gas and equi-luminous ground states consisting 
of 168 particles, d=2, x = 0.34523, 50 realizations. Clustering near the origin increases for increas- 
ing the constrained S{k) value, which are 1, 2, 3, and 4 respectively. The potential energy was 
minimized to within 10~^^ of its global minimum. 

FIG. 14: Ground-state configurations of 168 particles (a) S{k) = 2, x = 0.34523, and (b) S{k) = 
4, X = 0.34523. Particle clustering increases when constrained to higher S{k). 
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FIG. 15: Energy landscape and contour plot associated with minimizing Eq. ()Aip for the first wave 
vector. Three particles on a unit line, D = 0, x = 0.333. Ground state configurations are to the 
periodic arrangement. 



FIG. 16: Energy landscape and contour plot associated with minimizing Eq. ()Aip for the first wave 
vector. Three particles on a unit line, D = 1, x = 0.333. Ground state configurations are a set of 
configurations at the $ = ring around the integer lattice points. 
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FIG. 5: Batten, Stillinger, Torquato 
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FIG. 6: Batten, Stillinger, Torquato 
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FIG. 8: Batten, Stillinger, Torquato 
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FIG. 9: Batten, Stillinger, Torquato 
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FIG. 10: Batten, Stillinger, Torquato 
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FIG. 11: Batten, Stillinger, Torquato 
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FIG. 12: Batten, Stillinger, Torquato 
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(a) S(k) = 2, X = 0.35 
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FIG. 14: Batten, Stillinger, Torquato 
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FIG. 15: Batten, Stillinger, Torquato 
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FIG. 16: Batten, Stillinger, Torquato 
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